Quantum percolation in two-dimensional antiferromagnets 
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The interplay of geometric randomness and strong quantum fluctuations is an exciting topic in 
quantum many-body physics, leading to the emergence of novel quantum phases in strongly corre- 
lated electron systems. Recent investigations have focused on the case of homogeneous site and bond 
dilution in the quantum antiferromagnet on the square lattice, reporting a classical geometric per- 
colation transition between magnetic order and disorder. In this study we show how inhomogeneous 
bond dilution leads to percolative quantum phase transitions, which we have studied extensively by 
quantum Monte Carlo simulations. Quantum percolation introduces a new class of two-dimensional 
spin liquids, characterized by an infinite percolating network with vanishing antiferromagnetic order 
parameter. 
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Percolation occurs in a variety of contexts, rang- 
ing from blood vessel formation to clusters of atoms de- 
posited on substrate surfaces. A fundamental question is 
whether the classical picture of permeating networks jij 
apphes as well to strongly fluctuating quantum systems. 
In this respect, low-dimensional quantum magnets offer 
an ideal playground to probe the interplay between quan- 
tum fluctuations and geometric randomness both theo- 
retically and experimentally 0,0]. 

Two electronic spins coupled by an antiferromagnetic 
Heisenberg interaction, H = J Si -82 (J > 0), combine 
into a singlet state with purely quantum character. Such 
singlets are the fundamental building blocks of many 
novel phases in quantum antiferromagnets. For instance, 
the ground state of a Heisenberg ladder is a coherent su- 
perposition of singlet coverings involving mostly nearest- 
neighbor spins. This resonating valence bond (RVB) state 
is realized in a number of antiferromagnetically cor- 
related materials for which magnetometry and neutron 
scattering measurements indicate the absence of mag- 
netic order down to the lowest accessible temperatures 
0. In contrast, the 5 = 1/2 Heisenberg antiferromag- 
net on the square lattice has an antiferromagnetically 
long-range-ordered ground state However, as shown 
in Fig. n^A), a square lattice can be decomposed into 
dimers and ladders in such a way that there are no two 
adjacent dimers/ladders. This suggests that, if lattice 
disorder can isolate some of these substructures, it may 
lead to quantum disordered ground states. 

Recently, the effect of site dilution on the square- 
lattice S = 1/2 quantum Heisenberg antiferromag- 
net has been probed experimentally in the compound 
La2Cui_p(Zn,Mg)p04, in which magnetic Cu^+ ions are 
replaced randomly by non-magnetic Zn^+ or Mg^+ ions 
|3j. The pure system {p ~ 1) displays long-range anti- 
ferromagnetic order Upon doping with static non- 
magnetic impurities, the order parameter is reduced un- 
til it vanishes at a critical doping concentration pc- Both 
experimental data ^ and numerical simulations H H 



indicate a magnetic transition at the geometric percola- 
tion threshold p* = 0.40725 P| b eyond which the sys- 
tem breaks up into finite clusters. An analogous result 
holds for the case of bond dilution, where a fraction p 
of bonds is cut at random in the system: again the on- 
set of disorder is observed at the percolation threshold 
Pc = P* = 0.5 4]. In both cases this implies that the 
magnetic transition is controlled by the underlying geo- 
metric transition. A completely different picture of the 
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FIG. 1: {color online) (A) A square lattice quantum anti- 
ferromagnet composed of dimers (blue bonds, b) and ladders 
(yellow bonds, b'). (B) A segment of ladder in the percolating 
cluster leads to local RVB states (yellow ellipsis). The spins 
involved in the RVB state are very weakly correlated to the 
rest of the cluster, causing an effective fragmentation of the 
percolating cluster into non-percolating subclusters. 



effect of disorder might instead emerge when resonating 
valence bonds dominate the ground state of the system. 
In fact, the spins involved in a local RVB/singlet state are 
weakly correlated to the remainder of the system, so that, 
from the point of view of magnetic correlations, we can 
roughly consider these spins missing from the geometric 
cluster they belong to, as sketched in Fig.Q^B). This im- 
plies that the percolating cluster is effectively fragmented 
by the formation of RVB/singlet states, and therefore 
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the magnetic transition is hindered with respect to the 
geometric transition, acquiring in turn the nature of a 
quantum phase transition. Nonetheless, according to this 
argument the magnetic quantum transition appears as a 
renormahzed percolative transition, with an effective ge- 
ometry of the clusters dictated by quantum effects. This 
scenario is referred to as quantum percolation. 

In this work we investigate the occurrence of quan- 
tum percolation in the S — 1/2 quantum antiferromagnet 
on the square lattice with inhomogeneous bond dilution, 
modeled by the Hamiltonian 



— JbSl^b ■ S2,b + Jb'Si^b' ■ S2,b' 



(1) 



Here S are S — \/2 spin operators, the bonds b are dimer 
bonds in Fig. ^B), and the bonds h' are the inter-dimer 
bonds. Jb , Jb' > are random antiferromagnetic cou- 
plings drawn from a bimodal distribution taking the two 
values J ('on') and ('off'), with different probabilities 
of the 'on'-state, p{Jb — J) = P and p{Jb' — J) — P' ■ In 
the special case P — P' the homogeneous bond-dilution 
problem is recovered: the percolation threshold coincides 
with a magnetic order-disorder transition at the classical 
critical value Pc = P^ ^ 0.5 2]. In contrast, for P ^ P' 
one encounters inhomogeneous bond-dilution. In partic- 
ular, at P = 1 the system is an array of randomly coupled 
dimers, and in the opposite limit P' = 1 the system is a 
set of randomly coupled ladders. 
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FIG. 2: Phase diagram of the inhomogeneously bond-diluted 
S — 1/2 antiferromagnet on the square lattice. The colored 
area indicates the quantum-disordered region in which the 
system has developed an infinite percolating cluster, but the 
magnetization m of the system vanishes because of quantum 
fluctuations. The angle 6 parametrizes the critical curve. 

We determine the geometric percolation threshold for 
the inhomogeneous bond-dilution problem by generaliz- 
ing a very efficient algorithm recently developed for ho- 
mogeneous percolation . The resulting phase boundary 
is shown with filled diamond symbols in Fig. [3 A contin- 
uous transition line, with constant critical exponents of 
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FIG. 3: (A) Estimates of the quantum critical point as the 
value of the disorder showing the critical scaling (,x,^y ~ L at 
a temperature T = J/nL, plotted as a function of n. For n — 
8 the estimate clearly becomes temperature-independent. (B) 
Disorder- averaged square order parameter on the percolating 
cluster only as a function of the size Nc of the percolating 
cluster for the homogeneous percolation point {P = P' = 
0.5) and for a classical inhomogeneous percolation point {P = 
0.133, P' = 0.7). The temperature is T = l/(8Nc). The solid 
lines are cubic polynomial fits. 



2D percolation j^J, is found, connecting the percolation 
threshold for randomly coupled dimers, P = (P, P') = 
(1,0.382) with the percolation point P = (1/2,1/2) of 
the homogenously diluted systems, and terminating at 
the critical point for randomly coupled ladders, P = 
(0, 1). In this limiting case an infinitesimal concentra- 
tion of dimer bonds is sufficient to create a percolating 
cluster in the system. 

When classical spins {S = oo) are placed on each lat- 
tice site, long-range antiferromagnetic order follows triv- 
ially the onset of geometric percolation. In contrast, this 
picture changes drastically for = 1/2 quantum spins. 
Approaching the two limits of randomly coupled dimers 
(P ^ 1) and ladders (P' ^1), the formation of singlet 
states for the spins on the dimers (6-bonds) or on the lad- 
ders (5'-bonds) are favored statistically by geometry, such 
that quantum fluctuations are strongly enhanced. This 
contrasts with the case of homogeneous bond percolation, 
P — P', in which quantum effects do not change the clas- 
sical picture 4]. Therefore the crossover from homoge- 
neous to inhomogeneous percolation in the 5* = 1/2 an- 
tiferromagnet is accompanied by a continuous enhance- 
ment of quantum fluctuations. 

To determine the magnetic phase diagram in the quan- 
tum limit we use the Stochastic Series Expansion Quan- 
tum Monte Carlo (QMC) method based on the directed- 
loop algorithm [lOj. This technique enables us to ac- 
cess ultra-low temperatures on large Lx L lattices whose 
size is well inside the scaling range of the critical regime. 
To obtain finite-size scaling at T ^ we have carried 
out simulations on lattices as large as L = 64. For 
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strongly inhomogeneous systems, we have studied the 
T — > behavior by sitting at an inverse temperature 
/3J = 8L, at which the analysis of critical scaling leads 
to a temperature-independent estimate of the transition 
point, as shown in Fig. Ol^A). For systems closer to the 
homogeneous limit we found necessary to cool the sys- 
tem to even lower temperatures, which we can efficiently 
reach by a /3-doubling approach jj]. Finite temperature 
properties are studied on larger clusters with L = 200. 
Furthermore, simulations are also carried out at the per- 
colation threshold on percolating clusters of exact size 
grown freely from an initial seed (see Ref . 0] for a similar 
analysis). For each value of P and P' we average over 
100-^300 different realizations of the bond disorder (dis- 
order averages are indicated as (...)). The finite staggered 
magnetization {m?) = Y^3(5'(7r, tt)) / L"^ of an antiferro- 
magnetically ordered ground state is accompanied by a 
divergent spin-spin correlation length along the x and y 
direction, and ^y. The location of the magnetic tran- 
sition is obtained from the bond concentration at which 
the r = correlation lengths diverge linearly with lattice 
size, £,x,^y ^ L. 

The transition line for the quantum antiferromagnet 
and for the classical percolation problem are clearly dif- 
ferent. This shows that the magnetic phase transition 
in the case of inhomogeneous bond dilution turns into a 
quantum phase transition for sufficiently strong inhomo- 
geneities, leaving the percolation transition line at two 
multicritical points. In particular, between the two tran- 
sition lines there is an intermediate regime (colored area) 
in which geometrical percolation has been reached, such 
that there is an infinite percolating cluster Q . Nonethe- 
less, magnetic order is absent, as shown by the scaling 
of the order parameter in Fig. El^B). Here the disorder- 
averaged squared magnetization rric of the percolating 
cluster only is shown as a function of the cluster size. 
Unlike the case of the homogeneous percolation point 3 
P = P' = 0.5, the percolating cluster in the strongly 
inhomogeneous case shows a vanishing magnetization at 
the onset of geometric percolation, so that the addition of 
extra bonds is required for the magnetic order to set in. 
The quantum disordered regime of the phase diagram cor- 
responds to an infinite family of two-dimensional quan- 
tum spin liquids for which the antiferromagnetic Heisen- 
berg interaction does not lead to spontaneous symmetry 
breaking in the thermodynamic limit. 

Moreover, the phase diagram in Fig. |21 shows remark- 
able reentrant behavior. Decreasing P (P') from 1 to 
with fixed P' (P) in the range 0.45 < P' < 0.55 
(0.2 < P ^ 0.32), the system experiences first a magnetic 
transition from the spin liquid to the ordered antiferro- 
magnetic state. This is a clean example of the order-by- 
disorder mechanism, in which magnetic order is induced 
by randomness due to the geometric destruction of local 
RVB/singlet states. 

Unlike in other doped lattices do not ob- 



serve the occurrence of residual magnetic order due to 
the long-range effective interactions [l^ between spins 
not involved in a local RVB/singlet state. We ascribe 
this finding to a distinctive feature of bond dilution as 
opposed to other forms of disorder: in fact, when the 
cut of a bond does not lead to the local enhancement of 
quantum fluctuations, it leaves always two nearby spins 
not forming a singlet with their neighbors, so that each 
unpaired spin finds immediately its partner to build a 
longer-range singlet. 
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FIG. 4: (A) Scaling plot of (open symbols), (full sym- 
bols) for P = 1 (randomly coupled dimers). (B) Critical 
exponents z and v for the magnetic transition in the inhomo- 
geneously bond-diluted S=l/2 Heisenberg antiferromagnet. 
The angle 6 parametrizes the critical line in Fig. |5| Open 
(full) symbols refer to estimates obtained through the scaling 
of S,x (Cm)- The shaded areas correspond to the regions of 
parameters in which the transition is quantum. AH lines are 
guides to the eye. 

What are the critical exponents associated with this 
quantum phase transition? The finite-temperature quan- 
tum critical behavior of the correlation length ^ ~ T~^l^ 
determines the dynamical critical exponent z 0, [l3| . 
whereas the critical exponent v governs the divergence 
of the correlation length at the T = critical point, 
^ ~ |P — Pel Pc can be obtained as the value of 
disorder at which ^ exhibits a power-law behavior as a 
function of T, and it is in agreement with the estimates 
coming from the T = analysis described above. From 
the scaling relation 

i{-P\T) -}^F^{^T \P - P,r) (2) 
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the exponent v can be determined as the one reaUzing 
the best collapse of T^/''^{P-T) curves plotted versus 
TIP-PJ"^ for different P values, as shown in Fig. g^A). 
The values extracted for v and z are shown in Fig. EfB) 
as a function of the angle — tan~^[(l — P')/{1 — P)] 
which parametrizes the critical curve of Fig. [3 The v 
exponent is found to satisfy the Harris criterion for disor- 
dered systems v >2/ d = 1 within error bars. More- 
over, we observe that z and v are almost constant over a 
large portion of the critical line, taking the values z « 1.9 
and V « 4/3, and they both show a drop to z,v \ only 
near the edges of the critical line in the strongly quantum 
regimes (P ^ 1 and P' 1). Remarkably, such drop 
happens away from the multicritical points at which the 
percolation transition and the magnetic transition depart 
from each other. This implies that, for moderate quan- 
tum fluctuations, the magnetic quantum phase transition 
retains a percolative nature, in full agreement with the 
picture of quantum percolation, in which the main ef- 
fect of quantum fluctuations is a spatial renormalization 
of the percolating cluster. In particular, the correlation 
length exponent i/ is clearly consistent with the classical 
value z/ = 4/3 for two-dimensional percolation j2|] • The z 
exponent can also be related to the critical exponents of 
the classical phase transition. At and around the homo- 
geneous percolation point P = P' = 0.5, where the tran- 
sition is classical, the exponent z is the scaling dimension 
of the finite-size gap of the percolating cluster diverging 
in the thermodynamic limit, A ~ L^^. Since the num- 
ber of sites Nc of this cluster diverges as Nc ^ with 
D = 91/48 if the finite-size gap vanishes as A ~ A'^ 
it is easy to see that z = pD. For clusters with the ge- 
ometry of the chain or of the square lattice, p = 1 
Extending this result to the fractal percolating cluster we 
obtain z=£'=1.89.., in agreement with the QMC result. 

The concept of quantum percolation has previously 
been invoked in the context of weakly interacting elec- 
tronic systems, such as in metal-to-insulator transitions 
of granular metals 16], and in the integer quantum Hall 
effect of semiconductor heterostructures [ig. In these 
examples, quantum mechanics enters mostly through 
single-particle phenomena, such as localization effects of 
single-electron wave-functions and quantum tunneling. 
In contrast, the quantum percolation discussed in the 
present work occurs in a strongly correlated electron sys- 
tem. Its quantum aspects reside in the formation of local 
many-body quantum states, i.e. dimer singlet states or 
RVB states, and it ultimately has the nature of a col- 
lective quantum phase transition. The specific exam- 
ple of inhomogeneous bond dilution of the square-lattice 
Heisenberg antiferromagnet offers a remarkable realiza- 
tion of such scenario, with continuous tuning of the quan- 
tum nature of the percolative transition in the system. A 
continuous modulation of quantum effects at a percola- 
tion transition can also be achieved by progressively cou- 
pling two square-lattice layers which are site-diluted at 



exactly the same lattice sites [l8j. In that case, nonethe- 
less, the estimated critical exponents for the magnetic 
quantum phase transition differ substantially from those 
of classical percolation and are much closer to those of a 
quantum phase transition in a clean system, suggesting 
that quantum effects alter the percolative nature of the 
transition. 

From the experimental point of view, if different non- 
magnetic ions mediate the intra-dimcr interactions (6 
bonds of Fig. H^A)) with respect to the inter-dimer ones 
(6' bonds of Fig. H^A)) in a square lattice antiferromag- 
net, selective chemical substitution of these two types of 
ions offers a way to realize inhomogeneous bond disor- 
der. Extensions of the picture given in this paper would 
be necessary to mimic a realistic experimental situation. 
Nevertheless, the effects leading to a quantum correction 
of the percolation threshold should be robust, and they 
should persist in presence of more general models of in- 
homogeneous bond disorder. 
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